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Abstract 

The Standard Model fit prefers values of the Higgs Boson mass mn that are below 
the 114 GcV direct lower limit from LEP II. The discrepancy is acute if the 3.2a 
disagreement for sin^Oy^P^ from the two most precise measurements is attributed to 
underestimated systematic error. In that case the data suggests new physics to raise 
the predicted value of mn- One of the simplest possibilities is a Z' boson, which would 
generically increase the prediction for mn as a result of Z-Z' mixing. We explore 
the effect of Z-Z' mixing on the mn prediction, using both the full data set and the 
reduced data set that omits the hadronic asymmetry measurements of sin^^^^^, which 
are more likely than the leptonic asymmetry measurements to have underestimated 
systematic uncertainty. 



1. Introduction 

The Standard Model fit of tlie precision electroweak data, reviewed below, lias a less 
than robust confidence level, CL(x^, A^) = CL(17.2, 12) = 0.14, as a result of the enduring 
3.2cr discrepancy between the two most precise measurements of the effective leptonic weak 
mixing angle, sin^^^^^, from the polarization asymmetry A^ij and the front-back h quark 
asymmetry J^p^. Since the SM fit is relied on to provide guidance on the mass of the Higgs 
boson, it is relevant to consider the consistency of the sector of measurements that predict the 
value of vTiu- In this sector the problem is more severe, with CLlx^, N) = CL(14.1, 7) = 0.05. 
The discrepancy between Alr and Ap^ is reflected in a 3.2a discrepancy between the three 
leptonic asymmetry measurements, A^r, Ap^, y4^(P^), and the three hadronic asymmetry 
measurements, A'^^^, Ap^, and Qpb, and in the poor for the combination of all six 
asymmetries, CL(11.8,5) = 0.037. P These discrepancies could be statistical fluctuations, 
evidence of new physics, or the result of underestimated systematic uncertainty. If they are 
due to new physics, we cannot extract the Higgs boson mass mn from the precision data 
without flrst specifying the nature of the new physics. 

It might appear that the viability of the SM could be enhanced if the discrepancies are 
attributed to underestimated systematic uncertainty, in particular, in the hadronic asymme- 
try measurements, which share challenging, common experimental and theoretical systematic 
uncertainties. Indeed, if that is assumed and the three hadronic asymmetry measurements 
are omitted from the flt, the confidence level increases from 0.14 to 0.78, but a new problem 
emerges: the remaining measurements, dominated hj Apji , rriiy, and m^, predict rriH = 50 
GeV, with only a small probability, CL{mH > 114) = 0.03, for mn in the region mn > 114 
GeV allowed by the LEP II direct search limit. [2] Therefore this scenario also suggests new 
physics, in this case new physics to raise the predicted value of mn, and once again uih 
cannot be extracted from the data without specifying a model for the new physics. 

With this motivation several models of new physics have been considered to raise the 
the predicted value of rriH in the fit with hadronic asymmetries excluded, including light 
sneutrinos and gauginos,^ a fourth family of quarks and leptons,[l] and mixing with heavy 
vector-like leptons[5]. In this paper we consider mixing of the SM Z boson with a heavy 
Z' boson associated with a new Abelian symmetry with generator Qx, a simple extension of 
the SM that can raise the predicted value of mn- The mechanism is easy to understand: a 
heavy Higgs boson makes a negative contribution to the p parameter, p = my^/m|cos^6'vy, 
while mixing of Z with a heavier Z' shifts mz downward, causing p to increase so that 
the two effects tend to cancel. This possibility has been explored by Ferroglia, Lorca, and 
van der Bij[n] for the reduced data set with A^p^ excluded, for Z' bosons coupled to weak 
hypercharge Y and to B — L, the difference of baryon and lepton number. Our results agree 
qualitatively with theirs but differ in detail, both in the formulation of the Z' model and in 
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the implementation of the experimental constraints. In our approach but not in theirs the 
Z-Z' mass matrix is generated by Higgs bosons in the conventional wayj^ This theoretical 
difference has experimental consequences which are discussed below. They fit a truncated 
data set that captures the principal features but differs in detail from our fits, which are based 
on the complete EWWGpQ data set, use ZFITTER[7] to compute the radiative corrections, 
and include the largest experimental correlations as given by the EWWG. In addition, we 
impose the constraints on Z' production extracted by Carena et a/.|8j from the LEP II 
bounds for BSM contact interactions, [9] which we find are stronger than the precision EW 
constraints in parts of the Qx parameter space. We also impose the more recent constraints 
on Z' models obtained by the CDF collaboration, [TU] which are stronger than the LEP II 
bounds for some of the Qx parameter space if the Z' coupling constant is sufficiently smallj^ 
in particular, smaller than electroweak strength. Fits both with and without the hadronic 
asymmetries are presented. 

Following [m [12] we consider a class of models in which (1) the Z' receives its mass 
from a heavy SM singlet Higgs boson if', (2) the new gauge group U{\)x is required to 
be anomaly free with matter fields restricted to three SM generations augmented only by 
three right-handed neutrinos, and (3) the Qx charges of the SM fermions are independent 
of generation. It then follows [HI [12] that Qx must act on quarks and leptons like a linear 
combination of SM hypercharge Y and B — L, saj]^ 

Y B — L 

Qx = cos9x— + sinOx — - — • (1) 

The SM fermions and the W^, Zq bosons obtain their masses from the usual SM Higgs 
boson, which to preserve Qx gauge invariance must also be assigned Qx charge cos6'x|- + 



sinOx^— with its usual SM Y and (vanishing) B — L charges. In this approach, in order 
for there to be Z-Z' mass mixing, the SM Higgs boson must have nonvanishing Qx charge, 
QxH 7^ 0, requiring \9x\ 7^ ti/2. This contrasts with the model of [B] in which Z-Z' mixing 
can occur even if Qx acts on SM quanta purely like B — L. Following [12], we use the 
freedom to define the SM i?-hypercharge gauge boson and the new singlet Z' so that kinetic 
mixing vanishes at the electroweak scale and Z-Z' mixing is completely described by the 
mass matrix for the relevant energies near the TeV scale. We assume the SM singlet Higgs 
boson H' has a very large vacuum expectation value, v' 3> f , and that the new vector boson 
is much heavier than the Z, mz' 3> mz- 



^The authors of [6 exploit the fact that the Higgs mechanism is not necessary to ensure renormalizability 
in the case of Abehan gauge bosons. 

thank Bogdan Dobrescu for bringing the CDF bounds to my attention. 

^ A peculiar third solution obtained in [11] is incorrect, because those authors apparently failed to consider 
the S'C/(3)^ X U{\)x anomaly. 
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Qx then coincides with the SM generators cos6'x|- + sin6'x'^^^ in its action on SM 
matter quanta but not in its action on BSM quanta such as H' . In this framework with 
6'x = there can be a class of "y-sequential" Z' bosons with charges Y' which are identical 
to the SM hypercharge Y in their action on SM quanta but differ in their action on BSM 
quanta, [m, [12] These models were described as "unaesthetic" in [TT], although with a caveat 
that did not survive the journal's editorial process but is reproduced here: "We are humbly 
aware that aesthetic judgements are subjective and time-dependent. The cockroaches of 
Troy probably did not understand why the Greeks were making so much fuss."|^ This class 
of models has an appreciable effect on the allowed range of mj^^ in the EW fits, and especially 
for the fit of the reduced data set. 

The issues raised by the precision EW data will continue to be important in the era of 
the LHC Just as it has played an important role in the development of the SM, the precision 
EW data can also help us to understand the discoveries that will be made at the LHC. But 
our ability to use the precision EW data for this purpose will be severely limited if we cannot 
resolve the ambiguity created by the anomaly. 

In section 2 we review the SM fit and predictions for the Higgs boson mass. In section 3 
we describe the class of Z' models to be considered. In section 4 we summarize the relevant 
LEP II constraints on Z' bosons, taken from [8] and and the more recent constraints 
from CDF.[10j In section 5 we present constraints from the fits to the precision electroweak 
data together with the LEP II constraints. Concluding remarks are given in section 6. In 
an appendix we show that for a F-sequential Z' the effect oi Z — Z' mixing on the SM fit 
can be fully represented by "pseudo-oblique" corrections, which account for both vacuum 
polarization and vertex corrections. 

2. The Standard Model Fit 

In this section we review the SM fit to the precision electroweak data. We use the data 
set and methodology of the EWWG[I] with one exception: we do not include the W boson 
width in our fits, since with a 2.5% error it is not a precision measurement in the sense of the 
others, which are typically measured to 0(0.1%) or better, and in any case it has no impact 
on the prediction for the Higgs boson mass. We include the largest experimental correlations 
as given by the EWWG and use ZFITTER[3 to compute the radiative corrections, including 
the two loop contributions to sin^^^^^ amd mvi/-[l3] Like the EWWG we perform a fit to 
the data, scanning over m^, Aa[^^j(m^), asijnz), and mj^ j^leaving the latter two parameters 
unconstrained. The fits use the most recent Fermilab measurement of the top quark mass, |14j 



"^See the footnote on page 10 of the scanned preprint of [TT] posted at http://ccdb4fs.kek.jp/cgi- 



|bin/img_index?197706199 

°We have verified that the fit is not affected by scanning on mz because it is much more precisely measured 



than the other observables. 
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Experiment 


EWWG SM Fit 


Our SM Fit 


Pull 


Alr 


0.1513 (21) 


0.1480 


0.1480 


1.6 


A'- 

^FB 


0.01714 (95) 


0.01643 


0.01642 


0.76 




0.1465 (32) 


0.1480 


0.1480 


-0.45 


^FB 


0.0992 (16) 


0.1038 


0.1037 


-2.8 


^FB 


0.0707 (35) 


0.0742 


0.0741 


-1.0 


^wIQfb] 


0.2324 (12) 


0.2314 


0.2314 


0.83 


mw 


80.398 (25) 


80.377 


80.374 


0.95 


Tz 


2495.2 (23) 


2495.9 


2495.9 


-0.3 


Ri 


20.767 (25) 


20.743 


20.744 


1.0 


<^h 


41.540 (37) 


41.478 


41.477 


1.7 


Rb 


0.21629 (66) 


0.21581 


0.21586 


0.65 


Rc 


0.1721 (30) 


0.1722 


0.1722 


-0.04 


Ab 


0.923 (20) 


0.935 


0.935 


-0.6 


Ac 


0.670 (27) 


0.668 


0.668 


0.07 


rrit 


172.6 (1.4) 


172.8 


172.3 


0.24 


Aa5(m|) 


0.02758 (35) 


0.02767 


0.02768 


0.29 


asimz) 




0.1185 


0.1186 




rriH 




87 


85 





Table 1: SM fit compared with the EWWG fit.[l5] 

mt = 172.6 ± 1.4 GeV. 

In table 1 our fit is compared with the most recent EWWG fit [15], where it is clear 
the two are virtually indistinguishable: with Tw omitted both yield /N = 17.2/12 (with 
correlations contributing —1.4). The difference for rrit is an artifact of our fitting grid, 
Arrit = .21 GeV, which has been overtaken by the increasing experimental precision. The 
consistency for all other quantities, one part per mil or better, shows that the coarseness of 
the rrit grid has not affected the quality of the fit. For the Higgs mass our central value is 
85 GeV, indistinguishable from 87 GeV obtained by the EWWG. 

We can see from table 1 that the less than robust confidence level of the SM fit, 
CL(17.2, 12) = 0.14, is a consequence of the 3.2 a discrepancy between the leptonic and 
hadronic asymmetry measurements. A^p^ is the measurement with the largest pull, 2.82(T, 
corresponding to a a nominal Gaussian confidence level of 0.0048. The significance of such 
an outlyer can be estimated by the probability that one of twelve independent measurements 
will fluctuate to > 2.82cr, which is 1 — (1 — 0.0048)^^ = 0.06, enough by itself to account for 
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mn (GeV) 


90% CL 


CL{mH > 114) 




Alr 




34 


10- <mH < 108 


0.07 




mw 




52 


15 < < 135 


0.10 




A'fb 




480 


170 <mH< 1000+ 


0.99 






^A^ 


50 


19 < mj^ < 126 


0.07 


mw 






52 


10 < m/f < 140 


0.10 




ApB ^ 


^Qfb 


480 


180 <mH < 1000+ 


0.99 



Table 2: Predictions for uih from various restricted sets of mi^-sensitive observables. The 
value of ruH at the minimum is shown along with the symmetric 90% confidence interval 
and the likelihood for m^^ > 114 GeV. 10— and 1000+ denote intervals extending below 10 
or above 1000 GeV. 

the less than robust confidence level of the global fit. If wc consider only the observables that 
are sensitive to uih, omitting Rij^Rc-, A^-,^ Ac^cfh which are not, the confidence level falls 
by a factor 3 to CL(14.1, 7) = 0.05, and A^p^ is again the leading outlyer with a pull of 2.80(7 
and nominal likelihood 0.0051. The probability for such an outlyer is then 1 — (1 — 0.0051)^ = 
0.035, which matches nicely with the 0.05 confidence level of the corresponding fit. For 
both of these fits the central value for the Higgs boson mass is mu = 85 GeV, but since 
CL[mH > 114) = 0.26 there is no significant confiict with the LEP 11 114 GeV lower limit. 

It is instructive to consider why A^p^ is the outlyer in these fits rather than Alr. The 
explanation is that there is an "alliance" between the leptonic asymmetry measurements and 
the W boson mass against the hadronic asymmetry measurements. The former favor very 
low values for the Higgs boson mass, rriH ^ 50 GcV, as opposed to the hadronic asymmetry 
measurements, which predict uih an order of magnitude heavier, uih — 500 GeV. In table 
2 the central value from a fit to A^p^ alone is rriH = 480 GeV, while Alr and mw predict 
values an order of magnitude smaller, mn = 34 and 52 GeV respectively. We see in table 
2 and figure 1 that these differences are significant: the 90% symmetric confidence intervals 
for mn from the leptonic asymmetries and from the nonasymmetry measurements do not 
overlap the 90% interval from the hadronic asymmetries. (Note that the 95% confidence 
level upper limit is just the upper limit of the symmetric 90% confidence interval.) 

The interpretation of the precision data then depends critically on how we interpret 
the discrepancy between the hadronic and leptonic asymmetry data. If it is a statistical 
fluctuation then the prediction for mn from the fit to the full data set is applicable, with 
central value 85 GeV and 95% upper limit 158 GeV. Since the Ap^ and Alr measurements 
both represent many years of careful work, it is also certainly possible that the discrepancy 
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Combined Asymmetries 




10 20 50 100 200 500 



mH(GeV) 

Figure 1: distributions as a function of rriH from the combination of the three leptonic 
asymmetries A^ji, Ap^, A^(Pt-) (sohd hne); the three hadronic asymmetries A^p^, Ap^, and 
Qfb (dashed hne); and the three m//-sensitive, nonasymmetry measurements, mw,^z, and 
Ri (dot-dashed hne). The horizontal hues indicate the respective 90% symmetric confidence 
intervals. 

is a genuine reflection of new physics, for instance, in the Zbb vertex. Because Rb agrees well 
with the SM prediction, it is straightforward to show that this hypothesis requires a very 
large (~ 20%) new physics contribution to the right-handed Zbb coupling. Popular models 
of new physics cannot readily explain the data, but there is not a no-go theorem and some 
possibilities have been explored. [16] If it is a genuine manifestation of new physics, then the 
new physics must first be known in order to use the precision data to predict tuh- 

The third possible explanation of the discrepancy is underestimated systematic error. 
The three leptonic measurements are comparatively straightforward. They are free of com- 
plications from QCD and hadronization, involve three quite different techniques with no 
common systematic uncertainties, and have a sensible x^? with sin^^^'^''^ = 0.23113 (21) 
and CL{l.Q/2) = 0.44. In contrast, the three hadronic measurements share challenging 
experimental and theoretical systematic issues, including heavy fiavor tagging, large QCD 
corrections, and, especially, reliance on hadronic Monte Carlo simulations to merge the QCD 
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niH (GeV) 


90% CL 


CL{mH > 114) 


All data 


85 


47 <mH < 158 


0.26 


A'^ps © A'ps ®Qfb excluded 


50 


2A<mH < 105 


0.03 



Table 3: Predictions for mn from fits with and without the hadronic asymmetries, 
corrections with the experimental acceptance. They combine to give sin 

2^y^ff = 0.23222 (27) 

with CL(0.02/2) = 0.99. The surprisingly small results from an underlying 14 parameter 
heavy flavor fit, with an even more surprising x^, CL(53, 91) = 0.9995. These small values 
could result from overestimated systematic errors, but then the significance of the discrep- 
ancy is exacerbated and the fit CL decreases: e.g., using just statistical errors for the three 
hadronic measurements the of SM fit increases to 20.2/12 and the CL falls to 0.06. 
Another possible explanation of the small values is that they reflect incompletely under- 
stood correlations, which would again point to the possiblity of underestimated systematic 
error. A more detailed discussion is given in the talk cited in [2]. 

Only future experimental results can help us to choose among the three possible expla- 
nations. In this work we focus on the third possibility, not because we know it to be more 
likely but rather to understand the consequences. It might appear at first glance that the 
problem for the SM would be resolved if the three hadronic asymmetry measurements are 
assumed to have underestimated systematic errors and are omitted from the fit. The fit 
is then robust, with the p-value increasing from CL(17.2, 12) = 0.14 to CL(5.63,9) = 0.78, 
but the prediction for mn becomes problematic, with central value itlh = 50 GeV, with the 
95% CL upper limit at niH = 105 GeV, and with only 3% probability in the region allowed 
by the LEP II lower bound, CL{mH > 114) = 0.031, excluding the SM at ~ 97%CL. The 
rriH predictions from the two fits are summarized in table 3. 

One might think since the CL and CL{mH > 114) are independent probabilities that 
their product would be a measure of the likihood that the data agrees with the SM both as 
to the precision measurements and the Higgs mass prediction. However the product is not a 
fair estimator because it does not reflect the many ways that two independent probabilities 
can yield a product of a given value P1P2. A better estimator is the combined probability 
Pc that the product of two independent, uniformly distributed probabilities is less than or 
equal to the product P1P2, which it is easy to show is given by 

Pc{PlP2) = PlP2{l - l0g(PlP2)). (2) 

The fit to all data then yields Pc'(0.14 ■ 0.26) = 0.16, little changed from the x^ likelihood 
alone, while the reduced fit yields a somewhat smaller value, Pc(0.78 ■ 0.03) = 0.11. Clearly 
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the consistency of the SM with the data is not improved by removing the hadronic asymmetry 
measurements but rather the nature of the problem changes while remaining no less severe. 

Like the discrepancy between the hadronic and leptonic asymmetry measurements, the 
conflict of the reduced data set with the LEP II bound also has the canonical three possible 
generic explanations: new physics, systematic error or statistical fluctuation. We focus here 
on the possibility that it is an indicator for new physics, and consider below a class of 
Z' models that can maintain the quality of the fit for the reduced data set while raising 
the rriH prediction into the allowed region above 114 GeV. 
3. Z' Models 

We follow the framework described in [HI [12] and explored in detail in [12] . Restricting 
the fermionic content to the three SM generations augmented just by three right-handed 
neutrinos and assuming that the Z' couples universally to the three generations, a new 
\J{\)x gauge group is constrained to act on the three fermion generations like an arbitrary 
linear combination of the SM hypercharge Y and B — the difference of baryon and lepton 
number, as in equation (l).[IIl [12] Our study is restricted to the case of a very heavy 
Z' boson. Referring to the original, unmixed heavy gauge boson as Z^' , we assume that the 

mass is generated primarily by a heavy SM-singlet Higgs boson H' with a large vacuum 
expectation value, v' ^ v = 246 GeV, and that irizoi ^ m^. In order to preserve the 
U{l)x gauge invariance of the SM Yukawa interactions, Qx must also act on the SM Higgs 
boson H as indicated by equation (1). The interaction of Z^' with H then gives rise to mass 
mixing between Z^' and the SM boson Z", resulting in the mass eigenstates Z and Z' . In 
this framework Z^ — Z^' mixing only occurs if 7^ 0. 

The upper 2x2 corner of the 3x3 — B — Z^' mass matrix can be block diagonalized, 
yielding the massless photon eigenstate and the residual 2x2 Z^ — Z^' mass matrix, written 
compactly as 



In equation (3) m^o = gzv/2 is the usual unmixed Z^ boson mass, where gz = g/ cosOw, g 
is the SU{2)l gauge coupling constant, and 6w is the weak interaction mixing angle. The 
quantity r is the ratio of the U{l)x gauge coupling gz' to gz, 



and rhz' is the ratio of the Z' mass to the Z mass. 




(3) 



mz' 



mzoi 



> 1. 



(5) 



mz' = 



niz 
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Diagonalizing the mass matrix the leading correction to the Z boson mass is 

5m| = -r2cos2^x^ (6) 
and the Z — Z' mixing angle 9m, defined by 

Z = cos Om Z° + sin Om Z^' (7a) 

Z' = cos Om Z^' + sin 6m ijb) 

is 

e. = (8) 

Per equation (5), equations (6) and (8) are correct to leading order in l/m|,. 

The effect of the shift in the Z boson mass on the radiative corrections can be encoded [1^ 
as a contribution to the oblique parameter T, 



aT, = (9) 

so that 

^ r^cos^6'x , . 

aTx = 10 

The negative sign in equation (6), that occurs because the "levels repel" in two body mixing, 
implies a positive sign for Tx, which then causes the EW fit to prefer larger values of the 
Higgs boson mass. 

The second manifestation of Z — Z' mixing on the radiative corrections is the shift in 
the Zff couplings due to the admixture of Z^' in the Z mass eigenstate. Including the 
oblique corrections the interaction is 

^/ = 9z{l + ^) g'fj^f (11) 

where / represents a quark or lepton of chirality L or R and g'jr encodes the Zff coupling, 

g'f = 9f + reMq{ (12) 

Here (?/ is the SM Zff coupling 

9f = til - Q^^w (13) 

where and q-^ are the weak isospin and electric charge of fermion /. The quantity xw 
in equation (13) is the oblique-corrected square of the sin of the SM weak mixing angle, 

xw -xw = z aTx, (14) 

1 — 2xw 

9 



Xw = sin^ 9i 



and Qx is the Qx charge of fermion /, 

ql^ = cos Y + sin Ox — ^ — (15) 

where y^,b-^,l-^ are respectively the weak hypercharge, baryon number, and lepton number 
of fermion /. In keeping with the approximation rhz' ^ 1 we kept only the leading term in 
9m in equation (12). 

For a given choice of 6x the effect oi Z ~ Z' mixing on the EW fit is determined by a 
single parameter, which we choose to be Tx- The shift in the Zff coupling, equation (12) 
is determined by 

e = r9M, (16) 
which, using equations (8) and (10) is determined by Tx, 

aTx 



cos 9 



X 



(17) 



The fits presented in the section 5 are obtained by scanning over Tx in addition to the 
four SM scanning parameters, m^, Aa^^^{mz),as{'mz), and mH- The value of Tx determines 
the "effective Fermi constant" of the Z' boson, defined as 



Li 7,' 



4^2 m|, 



Defining Gz analogously. 



Gz = (19) 



4v/2m| 



which is equal at leading order to the Fermi constant, Gz = Gp, we have 

Gz' = 7=r" = = — (20) 
Gz m^z' cos^ 9x 

Since Gz' is constrained by the LEP II bounds, for a given value of 9x we obtain constraints 
on Tx both from the EW fits and from the LEP II bounds. 

Before proceeding to the EW fits of the Z' models we briefly mention an amusing feature 
of the y-sequential models. In general the EW corrections from Z — Z' mixing include both 
an oblique correction Tx from the shift in the Z boson mass|17j and non-oblique corrections 
from shifts in the Zff couplings, equation (12), due to the Z^' component of the Z eigenstate. 
But for the case of a F-sequential Z' boson, 9x = 0, we find that both the oblique and non- 
oblique corrections can be fully parameterized by correlated "pseudo-oblique" parameters, 
S' and T', defined by 

r = -Tx (21a) 
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S' = -4(1 - xw)Tx (216) 

where now 

The precision EW fit for the F-sequential Z' boson model can then be extracted from the 
usual obhque fit by considering the line S = 4(1 — xw)T in the S, T plane with T < 0. 

This parameterization of the model immediately reveals that the value oi mu cannot 
be increased toward the TeV scale and into the domain of dynamical symmetry breaking, 
which requires positive T and small or negative Since the original model with Tx > 
yields the same physics as the pseudo-oblique representation with T' < 0, because of the 
compensating effects of e and 5", we also see that one cannot attach an absolute significance 
to the sign of weak isospin breaking. 

The equivalence of the two representations is explained by the fact that for the Y- 
sequential model the apparently non-oblique correction to the Zff couplings induces a 
rescaling of the SM hypercharge coupling constant, which in turn contributes to W3-B kinetic 
mixing parameterized by S. This is not true for the other models we consider with 6x 7^ 0, 
since the term proportional to B ~ L cannot be absorbed into a renormalization of any SM 
interaction. A derivation is presented in the Appendix. 

4. Direct limits on Z' bosons from LEP II and CDF 

Carena et a/. |8j have used LEP II boundsf^] on contact interactions to extract limits on 
a variety of Z' bosons. Their results constrain the Z' effective Fermi constant, that is, the 
ratio of Z' mass to coupling strength, and in some cases they provide a stronger constraint 
than the precision EW data. For the interesting class of models with < < 7r/2 the 
CDF collaboration has obtained bounds [TO] which are stronger that the LEP II bounds if 
gz' is sufficiently small, gz' ^ fi'z/4. Both direct and EW constraints are presented in the 
results presented below. In this section we summarize the LEP II and CDF constraints for 
the Z' bosons considered in the EW fits presented in section 5 below. 

The class of Abelian charges considered here, defined in equation (1), is equivalent, in 
the notation of Carena et al. to the group U{l)q+xu-, characterized by the parameter x which 
ranges from —00 to +00 — see their table I. It is easy to see that the corresponding charge 
is 

T — 1 4- — T 

Qx = -^Y+^-^{B-L) (23) 



so that their x is related to our 6x by 



tan^x = ^. (24) 

X — 1 



^See for instance figures (12) and (13) of the second paper cited in [5]. 
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Defining gz' as the corresponding coupling constant, the relation between the coupling con- 
stants, determined by gz'Qx = gz'Qx, is 

gz' = ^^zV2x2 - lOx + 17. (25) 

With this dictionary we can translate the bounds obtained in ^ to the notation used 
here. We will see below that the most interesting region in 6x for the precision fits is the 
first quadrant, < 6x < tt/2, corresponding to the interval 4 > x > 1. Within this interval 
the 95% CL limit is (see figure 1 of 0) 

— > (2.62 + 1.18x)TeV (26) 

gz' 

Using the dictionary, equations (24) and (25), and equations (4), (5), and (10), equation 
(26) implies a bound on Tx, 

- (30.1 + 15.5 tan 0x)2' 

valid for < 6x < 7r/2. However, it should be noted that for 6x very near it/ 2 there must 
be a stronger bound, as can be seen by comparing the limits on U{l)q+xu and U{1)b-xL in 
figure 1 of [8]. At x = 1 both of these f/(l)'s become B — L but in the figure the latter is 
bounded more strongly than the former. A stronger bound must then exist on Z' bosons 
with charge Qx for Ox near 7r/2, that could be extracted from the two lepton differential 
cross sections, which are however not publicly available. In the following we restrict ourselves 
to the conservative bound, equation (27) 

Although the models with the greatest effect on mn lie in the first quadrant, < < 
7r/2, it is also interesting to consider the case Qx = T^r, since it occurs in attractive left-right 
extensions of the SM and also because it is typical of models in the second quadrant (or, 
equivalently, the fourth quadrant, since only the sign of gz' ■ Qx is physical). For Qx = T^ji 
we have 6x = — 7r/4, which corresponds to x — oo in the notation of [H]. The bound for this 
case is not discussed in [8], and we have extracted it directly from the LEP II constraint on 
the RR contact interaction quoted in j9]. In addition to T^fi we will sample the following 
choices from the first quadrant: 6x = 0, tt/6, tt/3, and ll7r/24, for which the LEP II bounds 
on Tx and Gz' are given in table 4. 

The case of 6x = llvr/24 is interesting because we will see in the next section that it 
has an appreciable effect on the EW fit even though it is very near 6x = 7i"/2 corresponding 
to Qx = B — L, for which there is no Z — Z' mixing and therefore no effect on the EW fit. 
However the surprisingly large effect that is found on the EW fit is severely constrained by 
the direct limit from equation (27) quoted in table 4. 
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Ox 


Tx 


Gz' 





0.14 


0.0011 


7r/6 


0.084 


0.00088 


7r/3 


0.039 


0.0012 


ll7r/24 


0.0059 


0.0027 


-7r/4 


0.30 


0.0047 



Table 4: 95% CL upper limits on Tx and Gz' obtained from LEP 11 bounds on contact 
interactions. 



Ox 


r 


mz' (TeV) 


Tx 




0.27 


0.83 


0.11 





0.13 


0.70 


0.039 




0.081 


0.61 


0.019 




0.20 


0.78 


0.051 


7r/6 


0.098 


0.64 


0.019 




0.059 


0.54 


0.0095 




0.20 


0.75 


0.018 


7r/3 


0.098 


0.60 


0.0070 




0.059 


0.45 


0.0046 




0.24 


0.69 


0.0022 


ll7r/24 


0.12 


0.50 


0.0010 




0.072 


0.40 


0.00058 



Table 5: 95% CL upper limits on Tx and mz' from CDF[10J for given values of r = gz'/gz- 
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The corresponding bounds on Tx from the CDF collaboration [TU] are given in table 
5, translated from the notation of reference [8], which is followed in reference [10], to the 
notation used here. For each 6x in table 4, except 6x = — 7r/4 for which no bound is given 
by CDF, we present the implied limit on mz' and Tx for given values of the Z' coupling 
strength, parameterized as the ratio to the SM Z boson coupling, r = gz'/dz- Comparing 
tables 4 and 5, we see that the CDF bounds are stronger than the LEP II bounds for r ^ 1/4, 
becoming increasingly stronger as r decreases. The LEP II bounds depend only on the ratio 
gz'/'mz', i-e., on the effective Fermi constant G^/, independent of the value of gz', because 
with mz' ^ mz they arise purely from the interference of the high energy tail of the Z 
boson amplitude with the low energy tail of the Z' amplitude. At Fermilab for sufficiently 
small gz'i which corresponds to smaller niz' for fixed Gz'-, the data begins to be sensitive 
to the direct production term, i.e., the square of the Z' amplitude, giving rise to increased 
sensitivity and a stronger constraint. 
5. Electroweak Fits in Z' Models 

In this section we present fits to the precision EW data for the class of Z' models 
discussed in section 3, focusing on the effect oi Z — Z' mixing on the value of the Higgs 
boson mass obtained from the fits. We will use two statistical methods that illuminate the 
physics in different ways, because they answer questions that are different in detail though 
they are clearly related. The first is the classical frequentist method, which is used by the 
EWWG[I] and was used in the discussion of the SM fit in section 2. In this method the 
question is 'Without imposing any a priori knowledge, direct or indirect, of the value of 
niHi how well does the model describe the precision data and what prediction does the best 
fit make for the likelihood of different values of mjif?' The second approach, followed for 
instance in the analysis of Z' models in [6], might be termed "Bayesian," in the sense that 
it imposes external knowledge about tuh as a "prior" constraint on the fit and therefore 
assesses the extent to which the fit of the model to the precision EW data is consistent with 
that prior. This approach then answers a different question which might be stated as follows: 
'If the value of mu were known to have a specific value or to lie within a certain range, how 
well does the model fit the precision data?' 

Both questions are valid and interesting, and it is useful to see what each tells us about 
the compatibility of various values of tuh with the different models. We will consider each 
in turn, combining the constraints from the fits with the direct LEP II bounds on Z' bosons 
in table 4. 
5a. Frequentist Fits 

We first present the frequentist fits for the models listed in table 4. These fits contain 
thank the authors of [8] for correspondence on this point. 
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Data Set 




CL 


mn (GeV) 


mj/(95%) 


Ch{mH > 114) 


A 


16.5/11 


0.12 


85 


153 


0.24 


B 


5.63/9 


0.78 


50 


105 


0.03 



Table 6: SM fits A and B. mif(95%) is the usual frequentist 95% upper limit on niH ■ 

one more free parameter than the SM fits, which may be chosen to be the effective Fermi 
constant Gz' or equivalently Tx (see equation (20)). In the frequentist approach to the SM 
fits, the global minimum is determined by freely varying /S.a^l^{rnz)-,Ois{'i^z)-, and 
ruH- The 95% CL upper limit on mn is then determined by minimizing as ruH is varied 
away from its value at the global minimum, until the local minimum (i.e., for the given 
value of mH) has increased by Ax^ = 2.71, corresponding to the upper boundary of the 
symmetric 90% confidence interval for one degree of freedom, CL(2.71, 1) = 0.90. 

To extend this approach to the Z' models, we vary Tx in addition to rritT iS.a^^{rnz)-, 
OisiiTiz), and uih, to obtain the global minimum, reducing the number of degrees of 
freedom by one relative to the SM fit, and then vary both Tx and mn about the minimum. 
The 90% contour in the mn — Tx plane is then defined by Ax^ — 4.61 corresponding to 
CL(4.61,2) = 0.90, and similarly the 95% contour is at Ax"^ = 5.99. The allowed regions 
are then further constrained by the 95% exclusion limits on Tx from the LEP II bounds on 
contact interactions, table 4, which are superimposed over the contours from the precision 
data. 

Wc consider two data sets. Data set B excludes the three hadronic asymmetry measure- 
ments. Set A contains all of the measurements in table 1 except the jet charge asymmetry, 
Qfb, which we omit for simplicity. (To compute the correction to Qfb we would have 
to convolute the mixing-induced shifts in the Zqq couplings with the qq partial rates, and 
then unfold the result to obtain the effective value of sin^^^*^^.) Qfb has very little impact 
on the fit CL or on rriH'- figure 2 shows that A^p^ completely dominates Ap^ and Qfb in 
the x^ distribution, since the combined distribution is practically indistinguishable from the 
distribution of A^p^ alone. 

The SM fits to data sets A and B are summarized in table 6. m/^(95%) is the frequentist 
95% CL upper limit, at x^ = Xmin + 2.71. Fit A has an acceptable prediction for rriH but 
a marginal confidence level, while fit B has a robust x^ CL but a failed prediction for itlh- 
The Z — Z' model fits to data sets A and B are shown in table 7 and figures 3-7. The effect 
of Z — Z' mixing on data set A is to push A^p^ further from the experimental value than 
the 2.82(7 deviation of the SM fit, so that the minima for set A coincide with the SM and 
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Figure 2: distributions as a function of rriH for the combination of the three hadronic asym- 
metry measurements (sohd hne) and for each individually: Ap^ (dashes), Ap^ (dashdot), 
and Qfb (dots). 

increases rapidly away from the SM minimum. As shown in table 7, the minimum for 
data set A is then at Tx = for all values of 6x, implying zero mixing, 6m = 0. The 
value and the central value of niH are then identical to the SM fit, while the confidence level 
decreases since there is one fewer degree of freedom, from CL(16.5,11) = 0.12 to CL(16.5,10) 
= 0.09. For set B the fits favor nonzero but small mixing for models with 6x in the first 
quadrant, with modest decreases in the minimum and modest increases in ttlh, while the 
fit likelihoods decrease slightly. For T^ji as for all models with 6x in the second quadrant, 
the minimum is at zero Z — Z' mixing, except for 6x very near vr where Qx — 

Although the changes in the minima are modest at best, there is a substantial effect 
on the allowed range of Higgs boson masses in the case of data set B and a smaller effect 
for data set A. This can be seen in the 90 and 95% contours shown in figures 3-7. For 
the F-sequential model, 9x = 0, in the case of data set B the 90% contour extends to 
rriH = 260 GeV, a factor 2.5 beyond the 105 GeV 95% upper limit of the SM fitfl The 

*^Note that quadrants I and III in dx are physically equivalent, as are quadrants II and IV, since the 
overall phase of Qx is not physical because it can be compensated by the phase of gz'- 

^The 90% contour of the Z — Z' tit should be compared to the symmetric 90% confidence interval of the 
SM fit, whose upper boundary defines the 95% upper limit. The extreme of the 95% contour corresponds to 
the 97.5% upper limit of the SM fit. 
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Data Set A 






Data Set B 




0x 




CL 


niH (GeV) 


Tx 




CL 


niH (GeV) 


Tx 





16.5/10 


0.09 


85 


0.0 


5.56/8 


0.70 


61 


0.012 


7r/6 


16.5/10 


0.09 


85 


0.0 


5.39/8 


0.72 


70 


0.019 


7r/3 


16.5/10 


0.09 


85 


0.0 


5.22/8 


0.73 


70 


0.014 


ll7r/24 


16.5/10 


0.09 


85 


0.0 


5.00/8 


0.76 


70 


0.006 


-7r/4 


16.5/10 


0.09 


85 


0.0 


5.63/8 


0.69 


50 


0.0 



Table 7: Frequent ist fits for Z' models. 

LEP II upper limit from table 4, Tx < 0.14, does not impinge on the contours from the 
EW fit. The maximum reach in ttih occurs at Tx = 0.10, corresponding to Gz' — 8 ■ 10~^; 
if the Z' coupling were of electroweak strength, this would imply a mass hierarchy of order 
tTiz' /mz — 30 or mz' — 3 TeV, within the range of the LHC For fit A the increase in mn 
is smaller, with the extreme of the 90% contour reaching 230 GeV, a factor 1.5 above the 
SM 95% upper limit at 153 GeV. Although the SM fit of data set B predict smaller values 
of rriH than set A, with Z — Z' mixing set B is consistent with larger values than set A. 

Figures 4 and 5 show that these features persist for 6x = 7r/6 and 7r/3. The LEP II 
constraint on Gz' begins to limit the allowed region for data set B, both because the LEP 
II bound becomes stronger and also because as 9x increases toward n/2 the ratio of Gz' to 
Tx increases like 1/ cos'^Ox, as seen in equation (20). Even though 6x = llvr/24 is very 
near 7r/2, Qx = B — L, for which there is no mixing and no effect on the fits, there is still 
a significant effect on the allowed region in mn for data set B, shown in figure 6. However, 
the allowed range is severely constricted by the direct limit on Gz', which, as discussed in 
section 4, is likely to be even stronger than is shown in the figure. For Qx = T^pt the effect of 
Z — Z' mixing on mfj is weaker, as can be seen from the more vertical slopes of the contour 
lines above 114 GeV in figure 7, but there is no additional constraint from the LEP II upper 
limit on Gz' which is Tx > 0.30. 

To estimate the confidence levels for these fits to lie within the LEP II allowed regions 
for rriH and Gz' we use a Bayesian likelihood method that was developed in the second 
paper cited in |2] to compute CL{mH > 114GeV) for the SM fits. In that approach two 
Bayesian priors were introduced to convert unnormalized likelihood functions into normalized 
probability distributions from which confidence intervals could be extracted. The first prior 
is that rriH lies between 10 and 3000 GeV. The precise value of the limits is not critical 
since there is negligible support above 1000 GeV or below 10 GeV. The second prior is that 
logniH is the appropriate measure, a natural assumption since the EW corrections depend 
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Frequentist Fit A: 0x=O. Qx=Y Frequentist Fit B: 0x=O. Qx=Y 




mH(GeV) mH(GeV) 

Figure 3: 90% and 95% CL contours in the Tx — mn plane for frequentist fits to data sets A 
and B for the F-sequential model, Ox — 0. The right axis indicates the corresponding values 
of Gz' — Gz' /Gz per equation (20). The diamond indicates m^j, Tx at the minimum for 
the Z' model. (Note that for set A for all Ox the minimum is at = and the diamond 
is hiding on the x-axis.) The ellipse and dot-dash horizontal line display the central value 
and 90, 95% symmetric confidence intervals of mn for the SM fit (elevated above Tx — ^ 
only for clarity). The horizontal dashed line is the 95% CL upper hmit on Gz' extracted 
from LEP II data. The vertical dashed line is the LEP II 95% lower limit on mn- 



Frequentist Fit A: Oy^=-n/Q 
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Frequentist Fit B: eyr=-n/Q 
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Figure 4: 90% and 95% CL frequentist contours for Z' model with Ox — t^/^, as in figure 3. 
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Frequentist Fit A: By^^-n/S Frequentist Fit B: 9y^=-n/?> 




mH(GeV) mH(GeV) 



Figure 5: 90% and 95% CL frequentist contours for Z' model with 6x = 7r/3, as in figure 3. 



logarithmically on ttlh- 

This procedure was shown to be reasonable (or at least no more foolish than the con- 
ventional procedure) by the fact that it provided confidence intervals for niH similar to those 
obtained from the Ax^ method, e.g., for data set B the result was CL{mH > 114) = 0.030 
from the Bayesian likelihood method versus 0.035 from Ax^ with the data of the time. We 
now find CL{mH > 114) = 0.17 from the SM fit to data set A compared to 0.24 from Ax^, 
and 0.018 compared to 0.031 for set B. There is no reason that the two methods should 
agree precisely. An important difference is that the Ax^ method compares only the best 
fits at different values of ttlh, while the Bayesian likelihood method samples the complete 
distribution of scanned parameters (mj, Aa^, as) at each value of mn- 

The same method can be applied to the two dimensional distributions in mn and Tx- 
The natural measure for mn is again logarithmic. Since Tx represents a first order pertur- 
bation of new physics on the leading order SM, the natural measure for Tx is linear. We 
normalize the likelihood functions in the intervals < Tx < 0.25 and 10 < rriH < 3000 GeV, 



where again the results are insensitive to to the precise choice of limits, The results are 
tabulated in table 8, which displays the confidence levels for rriH > 114 GeV, both with- 
out and with the LEP II constraint on Tx- In addition we tabulate mj|/(95%), which for 
the Z' models is defined as the largest value of rriH on the 90% contour that is consistent 
with the LEP II bound on Gz'- Unlike mj7(95%) for the SM fits, the values quoted for the 
Z' models cannot be interpreted as refiecting a 5% probability for rriH > mH{95%). Again 
the impact oi Z — Z' mixing is greater for data set B, with the probability of the LEP II 
allowed regions increasing by an order of magnitude relative to the SM value, e.g., from 
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Tx > is a boundary condition imposed by Z — Z' mixing. 
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Frequentist Fit B: e^=llTT/24 
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Figure 6: 90% and 95% CL frequentist contours for Z' model with 9x = ll7r/24, as in figure 
3. 
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Figure 7: 90% and 95% CL frequentist contours for Z' model with 9x — — 7r/4, Qx — T^r, 
as in figure 3. The LEP II upper limit, Tx < 0.30, is off the graph. 
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Data Set 


Model 




CL(mj/ > 114) 


and CL{Tx > Tlep ii) 


A 


SM 


153 


0.17 






ex = o 


230 


0.44 


0.44 




Tl/6 


220 


0.42 


0.42 




71/3 


214 


0.39 


0.38 




ll7r/24 


202 


0.34 


0.15 




-7r/4 


220 


0.43 


0.43 


B 


SM 


105 


0.018 






Ox = 


260 


0.29 


0.29 




7l/6 


252 


0.28 


0.24 




7r/3 


221 


0.23 


0.12 




ll7r/24 


158 


0.14 


0.01 




-7r/4 


188 


0.18 


0.18 



Table 8: Fits of data sets A and B. For the SM m//(95%) is the usual 95% upper limit 
obtained by the Ax^ method. For the Z' models m//(95%) is the maximum value of rriH on 
the 90% frequentist contours (figures 3-7) that is consistent with the LEP II direct limit 
on Tx- The confidence levels CL(m// > 114) and CL(Tx > Tlep ii) are computed with the 
Bayesian likelihood method described in the text. The entries in the last column combine 
both the TTiH and Tx direct limits from LEP II. 
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Data Set A 


Data Set B 


Model 


r 


mij(95%) CDF 


mj^(95%) LEP II 


mH{95%) CDF 


mH{95%) LEP II 




0.27 


230 




260 




^x = 


0.13 


225 


230 


190 


260 




0.081 


206 




139 






0.20 


220 




224 




7r/6 


0.098 


207 


220 


163 


252 




0.059 


194 




143 






0.20 


214 




172 




7r/3 


0.098 


193 


214 


141 


221 




0.059 


188 




134 






0.24 


188 




134 




ll7r/24 


0.12 


183 


202 


127 


158 




0.072 


181 




124 





Table 9: Effect of CDF bounds on the Higgs boson mass from frequentist fits of data sets 
A and B. As in table 8, m//(95%) is the maximum value of niH on the 90% frequentist 
contours (figures 3-7) that is consistent with the CDF direct limit on Tx for given values 
of r = Qz/qz'- The values of mij(95%) required by the LEP II bounds on T^, which are 
independent of r, are shown for comparison. 

GL{mH > 114) = 0.018 for the SM to 0.29 for the F-sequential model. 

In table 9 we show the effect of the CDF constraints from table 5 on the frequentist fits 
of data sets A and B. In particular for each value of r = gz'/oz we display mj/(95%), defined 
as in table 8, as the largest value of ttih on the 90% contours (figures 3-7) consistent 
with the corresponding upper limit on Tx from table 5. There is no CDF constraint for 
Qx = T,R. 
5b. Bayesian Fits 

In the frequentist fits presented above we scanned over as a free parameter, with no 
prior assumption except the exceedingly mild prior, 10 < uih < 3000 GeV, that was used 
only to obtain the confidence levels in table 8 for the regions in the rriH^Tx plane allowed 
by the direct LEP II limits on itlh and Gz'- The EW precision data alone determines the 
outcomes of those fits, which make predictions about the value of mn that can be tested for 
consistency with the direct LEP II lower bound on itlh. In this section we follow a different 
procedure: we suppose that the Higgs boson has been discovered at a specific mass which is 
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imposed as a prior constraint on the fits and ask how well the models describe the precision 
data for that value of itlh- This is the approach followed in [6]. We refer to this procedure 
as Bayesian because it assumes a prior value for mu- 

In table 10 we present results for uih = 114, 225, and 300 GeV. Since uih is fixed these 
fits have one more degree of freedom than the corresponding fits in section 5a. For each fit 
we present the minimum x^, the corresponding confidence level, the change in relative 
to the SM, and the value of Tx at the minimum. When Tx at the minimum exceeds 
the LEP II limit tabulated in table 4, we instead evaluate the fit with Tx set to the limit 
(marked by asterisks in table 10), so that the quoted is then the smallest value consistent 
with the LEP II limit. 

For data set A the Bayesian Z' fits at mu = 114 GeV do not improve on the SM fit, 
and the confidence levels are lower than the SM CL. For rriH = 225 and 300 GeV the Z' fits 
of set A have larger CL's than SM fit but they are still unacceptably low, ^ 0.03 and ^ 0.01 
respectively. For data set B the Z' models have a greater effect on the fits, and in all cases 
they improve on the SM. For mn = 225 and 300 GeV the confidence levels of the Z' fits 
are larger than the SM CL's by one and two orders of magnitude respectively, and the Ax^ 
values are highly significant. The Z' fit for 6x = ll7r/24 is severely constrained by the 
strong LEP II limit on Tx- The Z' models with the greatest effect on the fits are in the 
range ^ 6x ^ 7r/3, with the effect for 6x — vr/3 restricted by the LEP II limit on Tx for 
the larger values of ttih. For ^ 6*^ ^ 7r/6 the confidence levels are quite acceptable all the 
way up to rriH = 300 GeV. The large values of Ax^ in table 10 are unambiguous evidence 
of the effectiveness of the Z' model for set B with mn = 225 and 300 GeV. 

Contour plots for these Bayesian x^ fits are shown in figures 8-12. For the Z' models 
we exhibit the 90% and 95% contours with the LEP II limits on mn and Tx superimposed. 
The 90 and 95% confidence intervals for the corresponding SM fits are indicated by the tick 
marks on the horizontal dot-dashed line, elevated above the x-axis for visibility. With the 
Bayesian prescription, these intervals mark the value of rriH at which CL(x^/A^) = 0.10 or 
0.05, with = 12 for the SM fit to set A and A^ = 10 for B. Similarly the Z' contour plots 
are the 90 and 95% trajectories in the {rriH, Tx) plane with A^ = 11 and 9 for A and B 
respectively. Table 11 presents 95% upper limits on ttlh from these fits, defined for the SM 
as the upper limit of the 90% symmetric Bayesian confidence interval and for the Z' models 
as the largest value of tuh on the 90% contour that is consistent with the LEP II limit on 
Tx- For data set A the 95% upper limits of the Z' models are lower than for the SM, while 
for set B the limits increases relative to the SM, by a factor ~ 2 to nearly 400 GeV for the 
F-sequential boson. A qualitatively similar conclusion was reached by Ferroglia et a/., [6] 
who used the statistical method that we refer to here as Bayesian. 

It is interesting to reflect on the differences in the run confidence intervals for the fre- 
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.021 
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0.0059* 
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7.1 
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21.6/11 


0.068 


0.028 


3.4 


11.6/9 


0.12 


0.23 


8.9 


300 


SM 


31.8/12 




0.0015 




28.7/10 




0.0014 






^x = 


24.6/11 


0.062 


0.01 


7.2 


11.7/9 


0.11 


0.23 


17.0 




7r/6 


24.9/11 


0.054 


0.01 


6.9 


11.8/9 


0.084* 


0.22 


16.9 




7r/3 


25.5/11 


0.025 


0.008 


6.3 


14.7/9 


0.039* 


0.10 


14.0 




ll7r/24 


27.8/11 


0.0059* 


0.003 


4.0 


22.0/9 


0.0059* 


0.01 


6.7 




-7r/4 


24.8/11 


0.10 


0.01 


7.0 


14.9/9 


0.15 


0.09 


13.8 



Table 10: Bayesian fits of data sets A and B assuming fixed values of rriH at 114, 225, and 
300 GeV. is the chi-square minimum and is the number of degrees of freedom. Tx is 
the value at the minimum unless it exceeds the LEP II limit in table 4, in which case the 
fit is evaluated at the LEP II limit, denoted by an asterix. CL is the confidence level and 
Ax^ is the x^ difference between the Z' model and the SM fit. 
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Bayesian Fit A: 0x=O. Qx=Y 



Bayesian Fit B: 6^=0. Qx=Y 




Figure 8: 90% and 95% CL contours in the Tx — mn plane for Bayesian fits, as defined in the 
text, to data sets A and B for the Y-sequential model, 9x = 0. The right axis indicates the 
corresponding values of Gz' = Gz'/Gz per equation (20). The diamond indicates mH,Tx 
at the minimum for the Z' model. (Note that for set A for all 6x the minimum is at 
Tx = and the diamond is hiding on the x-axis.) The ellipse and dot-dash horizontal line 
display the central value and 90, 95% symmetric (Bayesian) confidence intervals of rriH for 
the SM fit (elevated above Tx = only for clarity). The horizontal dashed line is the 95% 
CL upper limit on Tx extracted from LEP II data. The vertical dashed line is the LEP II 
95% lower limit on tuh- 

quentist and Bayesian fits of data sets A and B. Consider first the SM fits. In the frequentist 
fits the 95% upper limit for itlh (the maximum of the 90% symmetric confidence interval) is 
153 GeV for set A and 105 GeV for set B, while for the Bayesian fits the pattern is reversed 
with 143 GeV for A and 183 GeV for B. The difference is due to the smaller confidence 
level of the fits to set A, e.g., CL(16.54, 12) = 0.13 for the frequentist fit to set A, compared 
to CL(5.63, 9) = 0.78 for B. The greater reach in mn of the Bayesian SM fit to set B is a 
consequence of the higher confidence level of the fit at the minimum, which allows for a 
greater excursion in niH, even though niH at the minimum is smaller for B than for A. 
In the frequentist fits the Ax^ method is used to compute the confidence intervals. In that 
method one computes the change in from the x^ minimum, without regard to what the 
value of x^ actually is at the minimum, so there is no penalty for the larger minimum of 
set A, and the 90% interval reaches to larger mn because of the influence of the hadronic 
asymmetries. In a sense the Ax^ method is Bayesian, since it assumes the flt at the 
minimum as a prior and then estimates the likelihood for deviations from the minimum. 

For data set A with uih near 114 GeV the x^ minima occur at very small values of Tx, 
as can be seen in table 10. The resulting flts have essentially the same x^ minima as the SM 
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Data Set A 


Data Set B 


Model 


rriH (95%) 


Tx 


mH{95%) 


Tx 


SM 


143 GeV 




183 GeV 




^x = 


127 


0.01 


390 


0.13 


7r/6 


128 


0.01 


368 


0.084* 


7r/3 


128 


0.007 


300 


0.039* 


ll7r/24 


128 


0.003 


220 


0.0059* 


-7r/4 


124 


0.01 


300 


0.13 



Table 11: 95% upper limits on rriH with the corresponding value of Tx, from the Bayesian 
fits. 771^(95%) is defined as the largest value of uih on the 90% Bayesian contours (figures 8 - 
12) consistent with the LEP II upper limit on Tx- An asterix indicates that ttih is evaluated 
for Tx at the LEP II upper limit from table 4. 



Bayesian Fit A: 6^=-n/6 Bayesian Fit B: 6x=tt/6 




mH(GeV) mH(GeV) 

Figure 9: 90% and 95% CL Bayesian contours for Z' model with 9x — tt/G, as in figure 8. 
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Bayesian Fit A: Q^=t\/2, Bayesian Fit B: 0x=V3 




10 20 50 HOO 200 500 ' 10 20 50 100 200 500 



mH(GeV) mH(GeV) 

Figure 10: 90% and 95% CL Bayesian contours for Z' model with Qx — t^I'^i as in figure 8. 



Bayesian Fit A: ex=llV24 
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Bayesian Fit B: Bf_=\\-r\/Z^ 
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Figure 11: 90% and 95% CL Bayesian contours for Z' model with Qx = ll7r/24, as in figure 
8. 
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Figure 12: 90% and 95% CL Bayesian contours for Z' model with Qx = -7r/4, Qx = Ts^j, 
as in figure 8. The LEP II upper hmit, Tx < 0.30, is off the graph. 

fit, and since they have one fewer degree of freedom the CL is lower than the SM, reaching 
0.10 at a smaller value of mu. The opposite is true of the Bayesian fits to set B, for which 
the minima in the Z' models are appreciably lower than in the SM, occuring at larger 
Tx , and the robust CL allows for larger values of mu before the probability falls to 0.10. 
In the frequentist fits of Z' models using the Ax^ method, mu also reaches larger values for 
set B than for set A, but in that case the effect is due entirely to the improvements in the 
fit at larger mu from Z — Z' mixing and not at all to the more robust SM fit of set B. 

The effect of the CDF bounds on the predictions of the Bayesian fits for the Higgs boson 
mass is shown in table 12. There is little effect on the fits to data set A since the bounds 
on Tx from the EW fits alone are already very strong. In the case of data set B the CDF 
bounds have more impact, especially for smaller gz'- 
6. Discussion 

We have explored the effect of a conservative class of Z' models on the Higgs mass pre- 
diction from the EW fits, considering both the possiblity that the discrepancy is a statistical 
fiuctuation and that it is the result of underestimated systematic uncertainty In the first 
case we fitted essentially all the precision EW data, data set A, while in the second we con- 
sidered the data without the three hadronic asymmetry measurements, data set B. The fits 
show that the range of allowed values for van can be significantly expanded into the allowed 
region above 114 GeV for data set B while retaining an acceptable fit to the precision data, 
but for data set A the possiblities arc more restricted. In particular, because of the marginal 
confidence level of the SM fit to data set A, the Bayesian fits of the Z' models allow even 
smaller domains for mH than the SM. 
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Data Set A 


Data Set B 


Model 


r 


mjj(95%) CDF 


mjf(95%) LEP II 


mH{95%) CDF 


mH{95%) LEP II 




0.27 


127 




384 




^x = 


0.13 


127 


127 


254 


390 




0.081 


127 




211 






0.20 


128 




302 




7r/6 


0.098 


128 


128 


220 


368 




0.059 


123 




195 






0.20 


128 




235 




7r/3 


0.098 


128 


128 


196 


300 




0.059 


128 




188 






0.24 


128 




190 




ll7r/24 


0.12 


125 


128 


180 


220 




0.072 


124 




176 





Table 12: Effect of CDF bounds on the Higgs boson mass from Bayesian fits of data sets 
A and B. As in table 11, mij(95%) is the maximum value of uih on the 90% Bayesian 
contours (figures 8 - 12) that is consistent with the CDF direct limit on Tx for given values 
of r = gz/gz'- The values of mij(95%) required by the LEP II bounds on Tx, which are 
independent of r, are shown for comparison. 

All Data Data Set B 



N 




10 20 50 100 200 500 1000 10 20 50 100 200 500 1000 

mH(GeV) mH(GeV) 

Figure 13: fits to the complete data set of Table 1 and to data set B with the obhque 
parameter T > 0. The sohd line is the distribution for the obhque fit, with the corre- 
sponding value of T shown in the dot-dashed hne which is read to the right axis. The dashed 
line is the distribution for the SM fit. 
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Frequentist Fit: All Data Frequentist Fit: Data Set B 




mH(GeV) mH(GeV) 

Figure 14: 90 (solid line) and 95% (dashed line) contour plots of frequentist fits with obhque 
parameter T > 0. The position of the minimum is indicated by the diamond. 



This is likely to be a generic feature of the response of the two data sets to models of 
new physics, because it is typically easier to construct models that raise the prediction for 
niu than it is to address the peculiarities of the J^p^ anomaly, as would be necessary to raise 
the marginal confidence level of the fit to data set A. For instance, obliquely mediated weak 
isospin breaking can raise the prediction for mu toward the TeV scale without impairing 
the quality of the fits but cannot improve them significantly. Figure 13 shows that fits with 
the oblique parameter T > 0, which generically represents weak isospin breaking mediated 
by vacuum polarization, can flatten the distribution for large values of m^^, for both the 
fit to the complete data set with 12 doj and for the fit to data set B with 9 doj. The 
confidence levels of these fits, ~ 0.13 for the complete set and ~ 0.7 for set B, are very near 
the CL's of the corresponding SM fits at their minima. 

Figures 14 and 15 display the 90 and 95% confidence level contour plots in the rriH^ T 
plane for frequentist and Bayesian fits. In a reversal of what we found for the Z' models, the 
minimum for the all-data set is at nonzero T with elevated uih, while for fit B it coincides 
with the SM fit with 7 = and mn — 50 GeV. However the position of the x^ minima are 
not very significant, since the minima are extremely shallow, as is evident in figure 13. 

Weak isospin breaking is the basis for the effect of Z — Z' mixing on the mn predictions 
presented here. In the Z' models the effect is hmited to niH ^ 300 GeV because the shifts 
in the Z f f couplings are proportional to gz'dM, which in turn is proportional to Tx, 

9z' Om = gz ■ (28) 

cos 9x 

The limit on Tx is reached when gz' Om grows so large that the Z f f couplings deviate too 
far from their SM values, causing x^ to increase. 
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Bayesian Fit: All Data Bayesian Fit: Data Set B 




mH(GeV) mH(GeV) 

Figure 15: 90 (solid line) and 95% (dashed line) contour plots of Bayesian fits with oblique 
parameter T > 0. The position of the minimum is indicated by the diamond. 



The 3.2cr discrepancy in the SM determination of the weak mixing angle by leptonic and 
hadronic asymmetry measurements sends an ambiguous message. It reduces the confidence 
level of the SM fit and raises questions about the SM prediction for rriH, which averages 
a "bimodal" distribution of measurements favoring 50 {Alr, my/) and 500 GeV {A^p^), as 
shown in figure 1. The significance and meaning of the discrepancy can only be clarified by 
future experimental data. If for instance evidence is found for a large deviation of the right- 
handed Zhh coupling from its SM value, it would confirm the A^p^ anomaly as a genuine signal 
of new physics. If the discrepancy results from unresolved theoretical and/or experimental 
systematic uncertainties in the very challenging hadronic asymmetry measurements, the low 
SM prediction for inn that results from the remaining measurements (data set B) strongly 
suggests new physics to raise the predicted value into the experimentally allowed region. 

The simple class of models studied here offers a paradigm for how we can use the LHC 
together with the precision EW data to understand the underlying physics. If a Z' boson 
is discovered at the LHC, it will be important to compare its properties as measured at the 
LHC with the constraints of the EW fits. It would first be essential to study the leptonic 
and hadronic couplings of the Z' to determine whether it is in fact in the class of Z' bosons 
considered here and to measure the parameter dx- If so, the measurement of the coupling 
constant gz' and mass m^/ would determine the effective Fermi constant Gz'-, which in turn 
specifies the oblique parameter Tx that determines the Z' fits. With the Z' parameters 
known, the EW fit would make a prediction for the Higgs boson mass which could be 
compared with direct measurement of mij at the LHC. It is then important to ascertain how 
well such a program could be carried out at the LHC, [19] both in its original incarnation 
and after possible luminosity and energy upgrades. 
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This example illustrates the important role that the precision EW data can continue to 
play in the future. At the end of the day, when new physics has been discovered and studied 
at the LHC, we will want to consider how it affects the EW fit. A consistent explanation 
of both the high energy data and the precision EW data would be a powerful confirmation 
of the theoretical picture, just as high energy data together with the precision data have 
confirmed the SM as the correct zero'th order model. If the model used to describe the high 
energy measurements is not consistent with the precision EW data, it could mean that the 
model is wrong or that there is other undiscovered new physics affecting the EW observables. 

To realize the potential of such a program, combining both the high energy measure- 
ments and the low energy precision data, it is important to resolve the ambiguity that the 
A'^PQ anomaly casts over the current data. Future high statistics studies at a high intensity 
Z factory like the proposed Giga-Z project [20] could determine if the anomaly is a statisti- 
cal fiuctuation and would allow further study of the experimental systematic uncertainties. 
Additional work on systematic uncertainties with a theoretical component would also be 
essential, for instance, the merging of the radiative corrections to Z — > 66 with the experi- 
mental acceptance, which gives rise to a systematic uncertainty that is now very difficult to 
quantify. 

It is also possible that the LHC could illuminate the issue. For instance, in the framework 
of the models discussed in this paper, the discovery of a F-sequential Z' boson together with 
a 300 GeV Higgs boson would be compatible with the Bayesian fit to data set B but not 
to data set A. In general, discoveries at the LHC could favor a model that is consistent 
with one data set but not the other. The EW fit of the compatible data set could then be 
compared with the direct observations at the LHC to further constrain and test the model. 
In principle, if the ambiguities can be resolved and the precision can be improved, the EW 
fit could even be used to probe for additional new physics before it is directly observed, just 
as the radiative corrections to the rho parameter [21] enabled the prediction of the top quark 
mass scale before the top quark was discovered. 

Appendix: "Pseudo-Oblique" Parameterization of the F-Sequential Model 

We will show, as discussed in section 3, that the effect of mixing with F-Sequential 
Z' bosons on the EW fits discussed here can be represented by oblique parameters S',T', 
given by 



r = -Tx 

S' = - xw)Tx 



{Al) 



where from equation (17) with 9 



'X = 







aTx = e 




{A3) 
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The precision EW fit for the y-sequential Z' boson model can then be extracted from the 
usual obhque fit by considering the fine S — 4(1 — xw)T in the 5", T plane. 

Using equations (12 - 16) with dx — 0, the Zjj interaction, equation (11), is 

>C/ = 5z (l + ^) imL - Q^xw + 6|)/. (A4) 

Wc will show that the interaction can also be represented by the obliquely corrected SM 
Lagrangian 

C^f = 9z{l + ^) imL - q'xM- 
where x'^r has the usual oblique correction in terms of S' and T", 

x'w -xw^ 1 -^2xw ("^ ~ ^^^^ ~ xw)T'^ (A6) 
and the W boson mass also gets the ususal correction. 

To obtain equations (A5) and (A6) we substitute y — 2{q — t^) in equation (A4). The 
result to order 0(e) is 

^f = 9z{l + ^-e) J Ml - Q^ixw - e{l - xw))]/- {A8) 



Matching the prefactors of equations (A5) and (A8) implies that aT' — aTx — 2e = —ocTx, 
which establishes equation (Al). S' is then fixed by the remaining condition for the equiva- 
lence of equations (A5) and (AS), x'-^ — xw — e(l — xw), which implies 

" f^' XH.(1 - x^)T'] ^ -aTx f '^T^^ - + 1 - X J iA9) 



1 - 2xw V 4 / V 1 - '^xw 

and using equation (Al) yields the result, equation (A2), for S'. 

The corrections to the W boson mass now provide a nontrivial test of the equivalence 
of the 5", T' oblique representation with the original Z' model. Prom the original model the 
correction is 

^ = aTx. (AlO) 

Tfl-^y 1 ^Xyj 

Substituting the expressions for S' and T' in terms of Tx from equations (Al) and (A2) into 
the generic expression for 5my/, equation (A7), the result is precisely equation (AlO). For 
all other observables we consider (see table 1) the oblique corrections enter via the weak 
mixing angle xw or, in the case of the Z width F^, via the prefactor 1 + aT. We are then 
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guaranteed that oblique fits with the constraint S' = 4(1 — Xw)T' are equivalent to the fits 
of the original F-sequential Z' model. We have also verified the equivalence numerically by 
fitting the data using both representations. 

This work was supported in part by the Director, Office of Science, Office of High Energy and 
Nuclear Physics, Division of High Energy Physics, of the U.S. Department of Energy under Contract 
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